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Abstract 

Evaluation of nucleon decay modes and branching ratios in a non-minimal supersymmetric 
SO (10) grand unified theory is presented. The non-minimal GUT considered is the super- 
symmetrised version of the 'realistic' SO(10) model originally proposed by Harvey Reiss, and 
Ramond, which is realistic in that it gives acceptable charged fermion and neutrino masses 
within the context of a phenomenological fit to the low energy standard model inputs. Despite 
a complicated Higgs sector, the SO(10) 10 Higgs superfield mass insertion is found to be the 
sole contribution to the tree level F-term governing nucleon decay. The resulting dimension 
5 operators that mediate nucleon decay give branching ratio predictions parameterised by a 
single parameter, the ratio of the Yukawa couplings of the 10 to the fermion generations. For 
parameter values corresponding to a lack of dominance of the third family self coupling, the 
dominant nucleon decay modes are p — ► K + + and n — > K° + £7„, as expected. Further, 
the charged muon decay modes are enhanced by two orders of magnitude over the standard 
minimal SUSY SU(5) predictions, thus predicting a distinct spectrum of 'visible' modes. These 
charged muon decay modes, along with p — ► 7r + + and n — ► 7r° + v^, which are moderately 
enhanced over the SUSY SU(5) prediction, suggest a distinguishing fingerprint of this partic- 
ular GUT model, and if nucleon decay is observed at Super-KAMIOKANDE the predicted 
branching ratio spectrum can be used to determine the validity of this 'realistic' SO(10) SUSY 
GUT model. 



1 Introduction 



Nucleoli decay is by definition a baryon number violating process, and within the context of the 
standard model (SM) of particle physics is forbidden []. Yet there is a strong motivation for as- 
suming that baryon number violation occurs; particularly the fact that there is no baryonic analog 
of the electromagnetic gauge invariance (which guarantees the conservation of electric charge), 
the presence of a baryonic asymmetry in the Universe ||, and the violation of baryon number 
conservation by black holes J3J. Allowing baryon number violation then suggests that the SM is 
only a low energy effective theory, and as such the stability of the nucleon is brought into question. 
This view was further reinforced when the adoption of Grand Unified Theories (GUTs) as a SM 
extension appeared to explain a large number of questions left unanswered by the SM ||. The GUT 
scheme was introduced to attempt to unify the SM interactions under a single simple gauge group. 
Imposition of such an underlying GUT structure then provided a new mechanism by which baryon 
number violation could occur, and nucleon decay induced || [7J. 

This nucleon decay mechanism is due to the fact that for conventional GUTs, quarks and leptons 
are placed in the same multiplets of the GUT gauge group. The coupling of these multiplets to 
either gauge or Higgs boson representations then gives interactions that couple quarks to leptons, 
and below the GUT scale, produce effective operators that induce nucleon decay. These tree level 
operators are four fermion dimension 6 operators || g] built from two fermion-fermion-boson vertices 
by means of a gauge or Higgs boson exchange. As the low energy limit of the internal boson 
propagator is -g? {Mq is the mass scale at which the GUT is spontaneously broken), the four fermion 
interaction reduces at low energy to an effective four fermion vertex scaled by two inverse powers 
of Mq. It is this class of effective vertex that would mediate nucleon decay in non-supersymmetric 
GUT models. 

In the archetypal GUT - minimal non-supersymmetric SU(5) - first proposed by Georgi and 
Glashow |J, the unification scale is Mq ~ 5 x 10 14 GeV [|T0| , |TT|j , and predicts the most dominant 
decay mode to be p — > 7r° + e + with a partial lifetime of r p ~ 4.5 x io 29±L7 yrs |IT| . This is to be 

1 Baryon number is not conserved in the SM, as violation occurs in weak interactions via instanton effects and the 
triangle anomaly, but the rate is suppressed and also involves violation of 3 units of baryon number for 3 standard 
model generations, thus making it irrelevant to nucleon decay H. 
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contrasted with the experimental lower bound obtained from IMB-3 Collaboration, of t p > 5.5 x 10 



ML' 



yrs [13]]. Clearly the minimal SU(5) model predicts proton decay at too rapid a rate, thereby 
ruling it out as a realistic GUT candidate. Nucleon decay channels and partial lifetime predictions 
have been calculated for a variety of GUT models Hl4| , including non-minimal SU(5) (which includes 
a 45 Higgs rep in an attempt to predict the fermion masses), minimal and non-minimal SO(10), and 
an Eq GUT model. Unfortunately, all these models tend to fail on the basis of a unification scale 
Mq ~ (2 — 7) x 10 14 GeV, which implies an overly rapid nucleon decay rate as well as a prediction 
for sin 2 9w that is inconsistent with the high precision LEP measurements [15 |. 

As conventional GUTs are essentially condemned by these failings, attention has turned to the 
super symmetric GUT models (SUSY GUTs). Imposing supersymmetry - a symmetry that relates 
bosons and fermions - has the effect of doubling the particle content below the GUT scale, which 
results in the slowing of the SM gauge coupling running, and consequentially predicts a consistent 
gauge coupling unification at a higher scale. Thus, a SUSY GUT model not only addresses the 
matter of the consistency of the sin 2 9w prediction, but it also predicts a unification scale that 
is typically two orders of magnitude larger than that of conventional GUTs. This increase in 
the unification scale induces a suppression factor of order 10~ 8 in the decay rates of four fermion 
dimension 6 operators generated by boson exchange, placing the dimension 6 mediated nucleon 
partial lifetime predictions well beyond the experimental lower bound. However, with the advent of 
Super-KAMIOKANDE, even the decay mediated by the dimension 6 operators may be observable. 

Yet the extension to a SUSY GUT model permits a new operator, capable of being the dom- 
inant contribution to nucleon decay. This operator is a dimension 5 fermion-fermion-sfermion- 
sfermion effective operator |L6|, [L7| constructed from either two fermion-sfermion-Higgsino vertices 
or a fermion-fermion-Higgs and a sfermion-sfermion-Higgs vertex by means of a heavy colour triplet 
Higgino or Higgs exchange below the GUT scale. Such an operator then evolves down to the SUSY 
breaking scale, at which point the sfermions are 'dressed' by gaugino exchange to give an effective 
four fermion vertex that mediates nucleon decay. As the low energy limit of the dimension 5 oper- 
ator is scaled by nucleon decay via this operator generally dominates over those mediated by 
the conventional dimension 6 operators^. 

2 It is assumed that R-parity is invoked so to rule out dangerous dimension 4 operators. 
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Investigations of nucleon decay in a number of SUSY GUT models have been carried out, 



beginning with the supersymmetrised version of minimal SU(5) [T8|, [TJ|. Unlike its non-SUSY cousin, 
this model predicts the dominant nucleon decay modes to be p — > K + + and n — > K° + z/ M , and 
as the unification scale is Mq ~ 2.5 x 10 16 GeV the partial lifetime prediction is r p ^ K + +p ~ 10 29±4 

yrs m n, m n n, n m, m, 1^, p n, n, m, m, m, 13, n (with m H3 set to m g ). 

This prediction does not disagree with the experimental bound of r p ^ K + +p > 10 32 yrs (obtained 
from the water Cerenkov detector of the KAMIOKANDE Collaboration |$9|, [L3|]) due mainly to 
the large uncertainty resulting from the value of the Higgs/Higgsino colour triplet mass. Likewise, 
predictions for non-minimal SUSY SU(5), minimal SUSY SO(10) result in a marginal degree of 
compatibility with the experimental lower bounds on the partial lifetimes of the various nucleon 



decay channels [[L4]]. This marginal consistency suggests that an improvement on the experimental 
lower bounds could lead to either a rejection of nucleon decay via SUSY GUT generated dimension 
5 operators, or an observation of nucleon decay. Yet the uncertainties in nucleon partial lifetime 
predictions preclude model discrimination by rate. In order to distinguish the underlying SUSY 
GUT structure, the relative decay rate predictions within a model should be determined, and then 
used to identify the SUSY GUT candidate, once nucleon decay has been observed. 

With this strategy in mind, this paper presents the branching ratios for nucleon decay in a 
particular 'realistic' SUSY GUT model. The model chosen is a supersymmetrised version [[HJ of 
the non-minimal SO(10) GUT proposed by Harvey, Reiss, and Ramond |^TJ, which was constructed 
primarily to reproduce a consistent phenomenological fit to the observed SM fermion masses and 
mixing angles. This realistic non-minimal SUSY SO(10) model, like its non-SUSY counterpart, 
can be viewed as a sophisticated phenomenological one, as it supports a rather expansive Higgs 
sector that is responsible for the required Yukawa coupling texture. It will be shown that analysis 
of the various nucleon decay channels mediated by the dimension 5 operators of this model results 
in branching ratio predictions depending on a single parameter, with the branching ratios for some 
observable modes enhanced by factors of order 100 over the minimal SUSY SU(5) predictions. This 
in turn suggests that if nucleon decay is observed at Super-KAMIOKANDE, the p — > K° + fi + , 
p — > 7T° + fi + , and n — > tt~ decay channel may play a significant role in identifying the structure 
of the underlying SUSY GUT. 

In this paper, section 2 presents the non-minimal SUSY SO(10) model to be used, section 3 
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examines in detail the low energy quark-level effective lagrangian, while section 4 discusses the 
effective lagrangian at the hadronic level and presents the branching ratio predictions. Finally, in 
section 5 a discussion of these predictions and the conclusions that can be drawn from them is 
given. 



2 The Non-Minimal SUSY SO (10) Model 

As mentioned, this analysis is based on the non-minimal SO(10) GUT model of Harvey, Reiss, 
and Ramond f4]J, which has been explicitly constructed to generate a mass spectrum (including 
mixing angles) of the SM fermions from the GUT. An advantage of the choice of SO (10) as the 
gauge group is that the lowest dimensional chiral representation that accommodates the observed 
SM fermions is the 16, which allows for the assignment of one family of SM fermions plus a right 
handed neutrino, and does not include any mirror fermions. This in turn places constraints on 
the possible Higgs sector representations, since the fermion masses transform under SO(10) as 
16 x 16 = (10 + 126 ), g + 120 4 . (where S and A refer to the symmetric and antisymmetric parts 
respectively), implying that the allowed Higgs reps that couple to fermions to form SO(10) invariant 
Yukawa terms are the 10, 120 , and the 126 . The 'realistic' model of Harvey et al. is then 
constructed from the representations in such a way that SO(10) GUT is broken directly to the 
SM gauge structure of SU(3)c x SU(2) x U(l), and the GUT scale texture of Yukawa couplings 
incorporates the up quark mass matrix ansatze of Fritzsch |S] and the down quark and charged 



lepton mass matrix ansatze of Georgi and Jarlskog [43] in such a way that the Oakes relation 



44f results. The cost of such a model is the expansion of the Higgs sector well beyond that of 
most minimal models. The particle content of this model is given in terms of a 45 that is the 
adjoint of vector bosons, three families of fermions (16 1? 16 2 , 16 3 ), and a scalar sector composed of 
a 54, a complex 10, and three families of 126 ( 126 ] , 126 ?, 126 ? ) - all of which are required for a 
viable spectrum of fermion masses. Note that the phenomenologically observed mass spectrum can 
be produced without requiring the presence of the 120 rep, which has a Yukawa coupling to the 
fermions that is antisymmetric in generation indices (as the SM fermions are expressed in terms of 
a single chirality, and the spin fields occur in a product that is symmetric in Lorentz indices). 
The extension [HD| of this model to that of a SUSY SO(10) model is straight forward, as the 
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gauge, fermion and conjugate Higgs fields are converted into vector and chiral superfields, giving a 
superfield content of: 



Vector superfields : 45 

Chiral superfields : 10, 16_i , I6 2 , 163, 54, HL , W.2 , 126? 



This SUSY SO(10) model, like its non-SUSY counterpart, is distinguished by its sophisticated 
Yukawa texture, composed of the 10, 16, and 126 chiral superfield reps. Specifically, the model is 
defined in terms of its superpotential, and for the purposes of nucleon decay, the relevant terms of 
the superpotential for this SO(10) model are 

W = (AI61 x 16 2 + 516 3 x 16 3 ) x HSi + (al6 x x 16 2 + 616 3 x 16 3 ) x 10 

+c(16 2 x 16 2 ) x H62 + d{m 2 x 16 3 ) x T26 3 + M G 10 x 10 (1) 

with the superpotential expressed in terms of the SO (10) representations, and A, B, a, b, c, d as the 
undetermined GUT scale Yukawa couplings. 

The beauty of this globally supersymmetric model is that as 10 x 10 D 1 and 126 x 126 p 1, the 
only SO(10) invariant F-term that contributes to nucleon decay below the spontaneously broken 
SO (10) GUT scale is given by Figure |I[ 




Figure 1: The only F-term supergraph that contributes to nucleon decay. 

The key point here is that this superfield diagram has a Higgs/Higgsino mass insertion that 
involves only the 10 (the 120 reps are absent!), which implies that only the GUT scale Yukawa 
couplings of the 16's to the 10 are of relevance to the predictions of nucleon decay (i.e. a and b in 
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equation ([I]) are the only relevant couplings). In terms of the particle diagrams, the only tree level 
diagrams of concern are 




Figure 2: The two particle diagrams generated by the F-term supergraph of Figure [T]. 

Here the first dimension 5 diagram exhibits the exchange of a Higgsino of GUT scale mass 
(Mq), and so for momentum below the GUT scale the Higgsino propagator reduces to a factor 
of . The second diagram in Figure ^| involves the exchange of a GUT scale Higgs scalar whose 
propagator reduces to -gr, but due to the weighting of the sfermion-sfermion-Higgs trilinear coupling 
by one power of Mq, the resulting diagram also contributes to the dimension 5 operator with weight 
■jj^. Thus the effective dimension 5 operator, valid between the SUSY GUT scale and the SUSY 
breaking scale (assumed to be of order the electroweak scale) is a combination of both diagrams, 
and is represented by the effective vertex in Figure 0. 

3 The Dimension 5 Operators 

In order to evaluate these dimension 5 operator contributions to nucleon decay, the superpotential 
must be re-expressed in terms of the superfields corresponding to the SM content. This may be done 
in a two step process, which first involves the re-expression of the superpotential in a compact SU(5) 
notation, followed by a decomposition of the SU(5) superfields into their SM components. Such a 
decomposition can be used as the F-term of Figure [TJ is the only dimension 5 contribution to nucleon 
decay, and it relies only on the Higgs 10 of SO(10) which has the SU(5) decomposition 10 — » 5 + 5. 
(Note that the SU(5) decomposition of the 16 is 16— ►10 + 5 + 1.) Thus the superpotential terms 
that contribute to nucleon decay can be written as 

Wsu i5 ) = V2 X fM^ ba H 2p - \e a ^xf M^ d xt H{ + M G H?H 2a (2) 
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with Wgu(5) being valid at the GUT scale. Here M u and M D are 3x3 matrices in generation space 
that express in a compact form the Yukawa coupling texture expressed in equation ([[]). Below the 
SO(10) scale, the heavy Higgs superfield can be integrated out to give an effective superpotential 
(that is appropriate below the GUT scale but above the SUSY breaking scale). This effective 
superpotential is 



W 



eff 

SU(5) 



V2 



AM, 



G 



(3) 



Here the Greek indices a, /3, ... are SU(5) indices, the family indices are (a, b, c, d), and the index 
i runs from 1 to 3. Also, the Lorentz structure is suppressed, so to focus on the generation and 
SU(5) structure. Restriction to the tree level diagrams relevant to nucleon decay (Figure D then 
implies the Yukawa texture matrices for this effective superpotential are of the form 



M 



a 





M 



D 



(4) 



As it is assumed that this SUSY SO (10) breaks straight to the minimally supersymmetric 
standard model (MSSM), the superpotential can then be further decomposed into the SM quark 
and lepton superfields. Typically the decomposition for one family of left-handed SU(5) matter 
superfields are 



I 

—v 



X 



ai3 



1 

72 



" 


u 3 


-u 2 


-u 1 


-d 1 


-U 3 





Ut 


-u 2 


-d 2 


u 2 


-U x 





— u 3 


-d 3 


u 1 


u 2 


u 3 





-L 


d 1 


d 2 


d" 


L 






(5) 



where U{, Di and Li are the charge conjugations of the right handed SU(2) singlet up, down, and 
charged lepton fields. Substitution of this decomposition into the superpotential (equation (|3|)), 
results in an effective superpotential relevant to nucleon decay that is expressed in chiral superfields 
associated with the SM. The form of the effective superpotential in question is 



W. 



SM 



Jjk 



eff 



4 -L a M u ab U bl U ck M°D dj + c —U ai MV b L b U ck M°D dj 



rD ; 



Jjk 



2M, 



G 



(6) 



-e ijk (ulM"dl - d a Mgul)(u k c M% - d k c M^ d )} 
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From this superpotential it is clear that as of a result of the SU(2) content, there are two 
classes of F-terms; the (LLLL)^ and the (RRRR)j? terms (here the notation of reference |T7| is 
used to emphasise the SU(2) weak content of the operators). However as the (RRRR)^ terms 
are antisymmetric in generation indices (a, b, c) - due to the Bose statistics of superfields in a 
superpotential - their composition is such that they must contain either a charm or a top SU(2) 
singlet superfield. This superfield generation remains, to a first approximation, unchanged on the 
dressing of the operator by gluino or bino exchange at the SUSY breaking scale (SU(2) gaugino 
exchange is forbidden for these singlet superfields), and so the low energy four fermion operator 
contains either a charm or a top quark. This implies that the (RRRR)^ term contribution to 
nucleon decay is suppressed, leaving only the (LLLL)^ terms. The effective lagrangian relevant to 
nucleon decay is then obtained from the (LLLL) p term of the superpotential by the usual method 
(Ant = \ J ( dijWj I ^=^^0 + h.c.) - Ei 1 1^ 11=0 witn $ representing a chiral superfield, and 
and ip the scalar and fermionic parts), which results in vertices composed of two particles and two 
sparticles. These vertices are then renormalisation group evolved down to the SUSY breaking scale 
(~ 0(My/)) at which point the dimension 5 operator is converted to a dimension 6 operator via 
gaugino and Higgsino exchanges. This dressing is schematically shown in Figure [| 



f „ f 




Figure 3: The particle diagram for the nucleon decay operators after being dressed by the wino or 
charged higgsino exchange. 

Of all the gaugino and higgsino exchanges associated with the dressing of the (LLLL)i? dimen- 
sion 5 operators, the dominant contribution comes from the charged wino. The gluino, neutral 
gaugino, and neutral Higgsino exchange contributions to the dressed operator are suppressed as 



S 



their exchange is approximately generation diagonal and their contribution is thereby suppressed 
due to small Yukawa couplings of the first and second generation fields present in the dimension 5 
operators |20|, |2l] . The charged-Higgsino exchanges are also suppressed, due to their Higgs strength 
Yukawa couplings to the first and and second generation fermions, thereby leaving the charged wino 
as the dominant contribution to the loop integral. Complete calculations of the loop dressing by 
chargino eigenstates have been performed by Sakai [f[9| , and their implications for nucleon decay 
rates have been explored for minimal and non- minimal versions of SU(5) f29j. Here however, the 
simplifying assumption of wino dominance of the decay amplitude is invoked. Preforming the loop 
integration, results in a triangle diagram factor, that although it depends on mass eigen-values and 
mixing angles of the sparticles in the loop, can be approximated (in the pure charged wino exchange 
limit) by [0 

d A k 1 1 1 



9l 



a 2 



z(27r) 



k 2 



m 



k 2 



m 



w 



mi 



2n m~. — m 2 7 

" a 



In 



m w 
m 2 -, 



(7) 



mt 



m l~ m W 



ln^- 



w d W ""W' 

and so becomes a multiplive factor of the dressed four fermion operator. 

From the superpotential (equation (0)) the effective lagrangian for the dressed quark level op- 
erators of Figure @) can be obtained, and with the use of equation (H) it has the form 



a 2 



R s R L M^ b M^e ijk {u l a d{){d k c v d ){f{u c , l d , m w ) + f(u a , d b , m w )} 



2nM G 

+ (d a ui)(u k c l d ){f(d c ,is d ,m w ) + f(d a ,u b ,m w )} 

+ (uldl)(u k J d ){f(u c ,d b ,m w ) + f(d a ,v d ,m w )} (8) 
+ (dy d )(dlu k c ){f(u a ,l d ,m w ) + f(u b ,d c ,m w )} ] + h.c. 

Here the Rs and Rl are the short and long range renormalisation factors. The short range renor- 
malisation accounts for the renormalisation effects from the SO (10) to the SUSY breaking scale, 
while the long range factor is from the SUSY braking scale to a low energy scale (assumed here to 
be 1 GeV). Rs can be shown to be generation independent, and can be taken to be |2y, [2TJ] 



Rf 





-4 




-3 




a 3 (m s ) 


9 


a 2 (m s ) 


2 












«G 



0.91 



(9) 
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where ms is the SUSY breaking scale (which here, has been set to the electroweak scale mw)- 
The long range renormalisation is predominantly a result of QCD interactions between the SUSY 
scale and 1 GeV, and encompasses the renormalisation of the Yukawa couplings and anomalous 
dimension corrections to the four fermion operators. Again, following reference |21| , 



Ri 



a 3 (lGeV) 
a 3 (m c ) 



a 3 [m c ) 



-18 
25 



a 3 (m z ) 



-18 
23 



0.22 



(10) 



a 3 (m b ) 

This effective lagrangian, as written, is for four fermion operators with the quarks and leptons 
expressed in their gauge interaction eigenstates; this however is easily remedied by rotating from a 
gauge interaction to a mass eigenstate basis. Due to the mismatch in the rotations of the charge 
|, — 3, and —1 fields, the operators incur additional generation mixing. The rotation matrices 
appropriate to this model are obtained from the diagonalisation of the mass matrices which are 
defined in terms of the Yukawa coupling texture (as specified in equation (fj)) and the vevs of the 
various Higgs reps. For these AIw = \ Dirac masses (Iw denotes weak isospin), it is convenient to 
consider the SO(10) vev contributions in terms of their SU(5) content. The contribution of SU(5) 
vevs to the quark and charged lepton masses is as follows |4l| : 



< . . . > ~ 5 gives a contribution to the charge | mass 

< . . . > ~ 5 gives an equal weight contribution to the charge — | and —1 masses 

< . . . > ~ 45 gives a contribution of relative weight 1 : — 3 to the charge — \ and —1 masses 



The SO(10) Higgs vevs structure is then decomposed as 



< 10 > 



< 126] > 



< 126, > 



< 126, > 



r (along 5) 4 
t (along 5) 
s (along 45) 
q( along 5) 



p(along 5) 



where p, q, r, s, t are taken as complex vevs. The assumption of complex vevs allows for the genera- 
tion of soft CP violation through the process of symmetry breaking. Yet it is assumed that soft CP 
violation is not the sole source of CP violation in the model. Hard CP violation is also permitted 
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due to the fact that, unlike the non-SUSY model of reference [ffl|, the Yukawa couplings of equation 
([1]) are taken to be complex . 



As the masses in the low energy effective SUSY theory arise from the Yukawa couplings of the 
quarks and charged leptons to a single light Higgs doublet of AI\y = ~, the mass matrices can be 
formulated in terms of the vev of this light Higgs. From the SU(5) decomposition of the SO(10) 



Higgs vevs, this light Higgs vev is a linear combination of the doublets in the 10, 126 i, 126 2 . and 126 3 
(in the ratio \r + p\ : t : s : q), and so the GUT scale couplings appearing in the mass matrices can 
be read off. Yet as it is the mass texture at the SUSY breaking scale that must be diagonalised, the 
entries in these Yukawa coupling texture matrices at the GUT scale must be evolved down to the 
SUSY scale via the renormalisation group equations, as done in Dimopoulos, Hall, and Raby jR| for 
'realistic' Yukawa matrices of this form. The quark and charged lepton Yukawa matrices specified 
at the SUSY breaking scale are then the mass matrices that are diagonalised. From equations ([!]) 
and (|12"D the GUT and SUSY scale mass matrix textures of the quarks and charged leptons are: 



U 



D 



GUT scale texture 
P G 
Pg Q G 

Q G V g 

R G e llpG 

Roe-*" S G 

R G 
Rg -35 g 
T G 



SUSY scale texture 
P 



U 








D 



P 



Rtr*e 


R 




K Q 
Q V 

Re** 
S 


R 

-35 






8 d 
T 



(12) 



with the assignments 



P = ap + At V = bp + Bt 
R = ar T = br S = cs Q = dq 



(13) 



and the subscript G indicating entries defined at the GUT scale. Here, the zero entries in the 
mass textures are the result of accidental discrete symmetries, which if broken, allow the generation 
of non-zero entries by means of the renormalisation group equations as the mass matrices are 
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renormalised down to lower energies. This is indeed the case for the entries 8 U and 5d which occur 
due to the violation of a discrete symmetry at the GUT scale. 

Although both the Yukawa couplings and the SO (10) Higgs vevs are complex, thereby permitting 



both hard and soft CP violation, the entries in the mass matrix textures, as given in (J12 ) , have 
been rendered explicitly real by means of quark and charged lepton field redefinitions. It is then 
these (SUSY scale) matrices, with 8 real parameters and one phase, that are diagonalised and the 
mass eigenvalues fitted to the low energy data, following Dimopoulos, Hall, and Raby f45fl . The 
diagonalisation proceeds by means of unitary and biunitary transformations of the form [/^iag = 
V U UV^, ZJ^iag = V^DV^\ and L^iag = ViLV^, and in following the assumptions of reference |45 



that V » Q ~ 5 U » P and T >> S ~ 5d >> R, the approximate mixing matrices are of the 
form 





c 2 


«2 


0" 


v u = 


-s 2 


c 2 













1 




Cl 


-si 


0" 


v d L = 


Sl 


Cl 













1 




c 5 


S5 


0" 


Vi = 


SB 


c 5 













1 



1 

c 3 

-S3 
1 

c 4 

-s 4 





S 3 
C 3 



S4 








(14) 



with Si = sin#j and q = cos^. The angles define in these rotation matrices can then be determined 
by fitting the mass eigenvalues to the low energy data. Using the low energy input data of EH], the 
resulting phenomenological fit specifies the angles as 



si ~ 0.196 s 2 ^0.05 s 3 ~ 0.046 
s 4 ~ 0.0066 s 5 ~ 0.070 cos cp ~ 0.41±g^§ 



(15) 



This phenomenological fit may need some revision in view of the subsequent and more precise low 
energy data (especially in light of the recent improvement to the bounds on the CKM matrix entry 



V c b [0) but it is expected that any revisions will have small effects on our results, and so we 
continue to use the original fit. 
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4 The Hadronic Lagrangian and Branching Ratio Predic- 
tions 



With the mass eigenstate rotations defined by equations and (|T6|) , the low energy effective 
lagrangian of equation (H) can be explicitly evaluated in terms of the dimension 6 four fermion 
quark level operators. In focusing on nucleon decay, these quark level (qqql) operators can be 
restricted by energy conservation, to have a quark composition of only the u, d, and s quarks. This 
in turn results in only five distinct operators, namely 

O q (dudv a ) = eytfdVjf/i/J O q (sudv a ) = e ljk (s l u j )(d k iy a ) O q (udsv a ) = e ijk (u i d j )(s k u a ) 
O q (duul a ) = ey fc (dV)(u fc O O q {suul a ) = t ljk (s^)(u k l a ) (16) 

Thus, the effective Lagrangian, expressed at the quark level, can then be written as £ nuc i eon = 
12 C (qqql) O q (qqql). Here the C(qqql)'s are the coefficients of the distinct quark level operators 
O q (qqql), and are determined from equation by summing the coefficients of the equivalent four 
fermion effective operators, modulo Fierz transformations. By classifying nucleon decay in terms of 
its various allowed channels, the effective lagrangian for nucleon decay can be written as 

C(n,p — > 7r + j?i) = C(dudv i )O q (dudi'i) (17) 

C(n,p^n + lf) = C(duuli)O q (duuli) 

C(n,p^ K° + 1+) = C(suuk)O q (suuli) 

C(n,p — > K + + j?i) = C(sudi>i)O q (sudi i i) + C(dusi>i)O q (dusi i i) 

Yet these effective lagrangian contributions are in terms of quark level operators, and so inappro- 
priate for hadronic decay rate calculations. Instead, they must be converted to effective lagrangian 
contributions at the hadronic level, thereby permitting evaluation of the nucleon decay rates, which 
although calculated at the hadronic level, are expressed in terms of the coefficients of the quark 
level four fermion effective operators specified by equations (^) and (18). This conversion may be 



preformed using the chiral lagrangian techniques developed in references [47], [23], and which 



express general hadronic level decay rates in terms of coefficients of generic four fermion quark level 
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operators. The results of these decay rate calculations, in the notation of p|j], are as follows: 
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Here mg = m^ = t«a = 1150 MeV is the mass to be associated with the virtual baryon exchange, 
m n = m p is the nucleon mass, and T> = 0.81 and T = 0.44 are numerical factors. 

From these decay rates, it is then very simple to construct branching ratios, which have the 
advantage over decay rates in that most of the as yet unspecified factors hidden in the quark level 
operators C(qqql) divide out, leaving the branching ratios parameterised by the ratio of the GUT 
scale Yukawa couplings of the complex 10. The numerical predictions for the branching ratios of 
the most dominant proton and neutron decay channels, for a large range of this parameter, |, are 
presented in Figures |] and |5] respectively. For this numerical evaluation, the branching ratios are 
defined as 

T(N ^x + y) 



Br(N x + y) 



T(N — > anything) 



(19) 
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Figure 4: The branching ratios of the most dominant proton decay channels. 



where N represents either the nucleon, and the decay rate for iV — > anything has been taken as the 
sum of all the relevant decay rates listed in ([19]). 



5 Conclusions 

With the results of the analysis of nucleon decay in this non-minimal SUSY SO(10) model presented 
in Figures f| and |5|, a number of important conclusions can be drawn. The first and most significant 
point is that this model gives one-parameter predictions for all the relevant nucleon decay branching 
ratios. Once nucleon decay is observed through any two channels, the ratio | is determined, and 
all the remaining partial lifetimes of the proton and the neutron then have a definite prediction. 
As with the SUSY SU(5) models, this model predicts that for a large region of | parameter space, 
p — > K + + and n — ► K° + are the most dominant proton and neutron decay modes. This 
prediction could only be altered by a strong suppression of the GUT scale Yukawa coupling a relative 
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Figure 5: The branching ratios of the most dominant neutron decay channels. 



to the third family self-coupling b, as shown by the prominence of the p — > K + + P T , p — > n + + v T , 
n — > K° + P T , and n — > tt° + v T decay modes for | < 3 x 10~ 2 . Another striking feature is that 
for | > 10~ 2 the branching ratio predictions are insensitive to the actual value of the parameter, 
thereby implying a degree of robustness to the predictions, regardless of the relative importance of 
the 10 of SO(10) in the assumed form of the GUT scale texture. 

However, it is the relative strengths of some of the individual branching ratios that serve to 
identify this model, and in particular, it is the nucleon decay channels involving the ji + and the 
Pfi that are the distinctive fingerprints of this model. For both the proton and the neutron, the 
branching ratio predictions for channels involving the charged muon show a marked enhancement 
over corresponding predictions of minimal SUSY SU(5). Specifically, the branching ratio predictions 
for the p — > K° + fi + , p — > 7r° + fi + , and n — > 7r~ + relative to the dominant proton and neutron 
decay channels are enhanced over the minimal SUSY SU(5) predictions by factors of 50-500, 10-100, 
and 20-200 respectively (the ranges given in these enhancement factors are due to the uncertainty of 
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the minimal SUSY SU(5) predictions as quoted by |29], [25], [34]]). To a lesser extent, the p — > 7r + + u l 



and n — > ir + z/^ decay channels show a similar enhancement, but only by factors of 3.6 and 2.6 
respectively. Thus, these enhancements in the decay rate predictions result in branching ratios 
for this non-minimal SUSY SO(10) model that are both qualitative and quantitatively different 
from that of the SUSY SU(5) nucleon decay spectrum, thereby making this 'realistic' non-minimal 
model a testable candidate for a SUSY GUT extension to the standard model. The issue of testing 
the predictions of this model could be addressed at Super-KAMIOKANDE, provided that Super- 
KAMIOKANDE in fact observes nucleon decay. 

In sum, the distinctive tests of this realistic supersymmetric SO(10) GUT which arise from the 
consideration of nucleon decay come not from the actual decay rates or partial lifetimes of the 
nucleon, as the nature of the Higgs sector and the uncertainty of the Higgs and Higgsino colour 
triplet masses make the SUSY dimension 5 operator decay rate predictions uncertain. Rather, they 
come from the calculation of nucleon decay branching ratios. The fact that this realistic model 

predicts ratios of branching ratios f^fr^j , bXZk++v}) » and Br(n^Ar°+C) of order 20% ' shows 
the relevance of 'observable' channels such as p — > K° + p — > tt° + and n — > tt~ + fi + to the 
testing of models of GUT unification. (For related considerations involving mass textures induced 
by higher dimensional operators see P5| .) These enhanced branching ratio predictions are instead 
simply a result of the composition of the Higgs superfield sector, which is such that the GUT scale 
Yukawa couplings relevant to nucleon decay are not the full set of couplings that contribute to SM 
fermion mass generation. 

The results presented here may be seen as some of the possible implications of a viable SUSY 
GUT model, and any observation of p — > K°+fi + , p — > 7r°+ / u + orn^ tt~+[i + at a level significantly 
enhanced above the expected SUSY SU(5) predictions is an indication that the underlying structure 
of a realistic extension to the standard model is best described in terms of a SUSY GUT model with 
a non-minimal Higgs sector. Unfortunately, because only a partial set of the GUT scale Yukawa 
couplings is directly involved in the analysis of nucleon decay, whereas the light Higgs is a linear 
combination of contributions from the various SO (10) Higgs reps, the actual values of the GUT 
scale couplings remain undetermined and the texture unexplained, at least in this model. 
Acknowledgements 

The author wishs to thank B.A. Campbell for suggesting the problem and his subsequent supervision 



17 



and encouragement, as well as N. Rodning for the useful discussions that occurred over the course 
of the work. 

References 

[1] G. t'Hooft, Phys. Rev. Lett. 37 (1976) 8. 

[2] T.D. Lee and C.N. Yang, Phys. Rev. 98 (1955) 101. 

[3] A.D. Sakharov, JETP Letters 5 (1967) 24. 

[4] S.W. Hawking, Nature 248 (1974) 30. 

[5] For a review see P. Langacker, Phys. Rep. C72 (1981) 185. 

[6] H. Gcorgi and L. Glashow, Phys. Rev. Lett. 32 (1974) 438. 

[7] J.C. Pati and A. Salam, Phys. Rev. D10 (1974) 275. 

[8] S. Weinberg, Phys. Rev. Lett. 43 (1979) 1566. 

[9] F. Wilczek and A. Zee, Phys. Rev. Lett. B43 (1979) 1571. 

[10] H.Georgi, H.R. Quinn, and S. Weinberg, Phys. Rev. Lett. 33 (1974) 451. 

[11] M. Goldhaber and W.J. Marciano, Comm. Nucl. Part. Phys. 16 (1986) 23. 

[12] W. Gajewski, et. al., Phys. Rev. D42 (1990) 2974. 

[13] For a review of present limits see R. Barloutaud, Proceedings of the International Workshop 
on Theoretical and Phenomenological Aspects of Underground Physics (TAUP 91), edited by 
A. Morales, J. Morales, J.A. Villar (North-Holland, 1992) p. 522. 

[14] For discussions on various non-SUSY and SUSY model predictions, refer to G.G. Ross, Grand 
Unified Theories, Frontiers in Physics Series (Addison- Wesley, 1985), and R.N. Mohapatra 
Unification and Supersymmetry (Springer Verlag, 1986). 

18 



[15] P.Langacker and N. Polonsky, Phys. Rev. D47 (1993) 4028. 

[16] N. Sakai and T. Yanagida, Nucl. Phys. B197 (1982) 553. 

[17] S. Weinberg, Phys. Rev.D26 (1982) 287. 

[18] S. Dimopoulos and H. Georgi, Nucl. Phys B193 (1981) 150. 

[19] N. Sakai, Z. Phys. Cll (1982) 153. 

[20] S. Dimopoulos, S. Raby, and F. Wilczek, Phys. Lett. B112 (1982) 133. 

[21] J. Ellis, D.V. Nanopoulos, and S. Rudaz, Nucl. Phys. B202 (1983) 43. 

[22] W. Lucha, Nucl. Phys. B221 (1983) 300. 

[23] V.M. Belyaev and M.I. Vysotskii, Phys. Lett. B127 (1983) 215. 

[24] S. Chadha and M. Daniel, Nucl. Phys. B229 (1983) 105. 

[25] S.J. Brodsky, J. Ellis, J.S. Hagelin, and C. Sachrajda, Nucl. Phys. B238 (1984) 561. 

[26] N. Sakai, Nucl. Phys. B238 (1984) 317. 

[27] S. Chadha and M. Daniel, Phys. Lett. B137 (1984) 374. 

[28] J. Milutinovic, P. B. Pal, and G. Senjanovic, Phys. Lett. B140 (1984) 324. 

[29] B.A. Campbell, J. Ellis, and D.V. Nanopoulos, Phys. Lett. B141 (1984) 229. 

[30] L.E. Ibanez and C. Munoz, Nucl. Phys. B245 (1984) 425. 

[31] S. Chadha , CD. Coughlan , M. Daniel, and G.G. Ross, Phys. Lett. B149 (1984) 477. 

[32] P. Nath, A.H. Chamseddine, and R. Arnowitt, Phys. Rev. D32 (1985) 2348. 

[33] P. Nath and R. Arnowitt, Phys. Rev. D38 (1988) 1479. 



19 



[34] R. Arnowitt and P. Nath, Phys. Rev. Lett. 69 (1992) 725. 

[35] J. Hisano, H. Murayama, and T. Yanagida, Nucl. Phys. B402 (1993) 46. 

[36] J.L. Lopez, D.V. Nanopoulos, and H. Pois, Phys. Rev. D47 (1993) 2468. 

[37] K. Daum, Z. Phys. C62 (1994) 383. 

[38] T. Nihei and J. Arafune, Prog. Theor. Phys. 93, 665 (1995). 

[39] KAMIOKANDE-II Collaboration (K.S. Hirata et. al.), Phys. Lett. B220 (1989) 308. 

[40] B.A. Campbell, S. Davidson, and K.A. Olive, Nucl. Phys. B399 (1993) 111. 

[41] J. A. Harvey, D.B. Reiss, and P Ramond, Nucl. Phys B199 (1982) 223. 

[42] H. Fritzsch, Phys. Lett. B70, 437 (1977). 

[43] H Georgi, and C Jarlskog, Phys. Lett. B86 (1979) 97. 

[44] R.J. Oakes, Phys. Rev. D26 (1982) 1128. 

[45] S. Dimopoulos, L.J. Hall, and S. Raby, Phys. Rev. D45 (1992) 4192. 

[46] CLEO Collaboration, B. Barish et al, Phys. Rev. D51 (1995) 1014; 

M. Shifman, N.G. Uraltsev, and A. Vainshtein, Phys. Rev D51 (1995) 2217. 

[47] M. Claudson, M.B. Wise, and L.J. Hall, Nucl. Phys. B195 (1982) 297. 

[48] H. Murayama, and D.B. Kaplan, Phys. Lett. B336 (1994) 221. 



20 



